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Abstract 

We present an effective linear response approach to pump-probe femtosecond coherence spectroscopy in the well 
separated pulse limit. The treatment presented here is based on a displaced and squeezed state representation 
for the non-stationary states induced by an ultrashort pump laser pulse or a chemical reaction. The subsequent 
response of the system to a delayed probe pulse is modeled using closed form non-stationary linear response 
functions, valid for a multimode vibronically coupled system at arbitrary temperature. When pump-probe sig- 
nals are simulated using the linear response functions, with the mean nuclear positions and momenta obtained 
from a rigorous moment analysis of the pump induced (doorway) state, the signals are found to be in excellent 
agreement with the conventional third order response approach. The key advantages offered by the moment 
analysis based linear response approach include a clear physical interpretation of the amplitude and phase of 
oscillatory pump-probe signals, a dramatic improvement in computation times, a direct connection between 
pump-probe signals and equilibrium absorption and dispersion lineshapes, and the ability to incorporate co- 
herence such as those created by rapid non-radiative surface crossing. We demonstrate these aspects using 
numerical simulations, and also apply the present approach to the interpretation of experimental amplitude 
and phase measurements on reactive and non-reactive samples of the heme protein Myoglobin. The role played 
by inhomogeneous broadening in the observed amplitude and phase profiles is discussed in detail. We also 
investigate overtone signals in the context of reaction driven coherent motion. 



I. INTRODUCTION 

Femtosecond coherence spectroscopy (FCS) is an ultrafast pump-probe technique that allows the experimen- 
talist to create and probe coherent vibrational motions and ultrafast chemical reactions in real time pl-pof. In 
a typical pump-probe experiment, an ultrashort pump laser pulse is used to excite the sample of interest. The 
subsequent non-stationary response of the medium is monitored by an optically delayed probe pulse. Owing to 
the large spectral bandwidth available in a short laser pulse, it can generate non-stationary vibrational states 
in a molecular system as shown in Fig. 1(a). The subsequent nuclear dynamics modulates the optical response 
as detected by the probe pulse. Coherent vibrational motion in the ground state has been observed in crys- 
talline and liquid phase systems [p[p^, and in biological specimens having short excited state lifetimes such 
as bacteriorhodopsin HJ^, and myoglobin |lj]. For molecules that have long-lived excited states, however, the 
excited state coherence is dominant and has been identified in several dye molecules in solution ||l|,^,^ , in small 
molecules in the gas phase |^,0,|l^ , and in photosynthetic reaction centers . 

Apart from "field driven" coherence directly prepared by the laser fields, vibrational coherence can also be 
driven by rapid non-radiative processes. For example, if we consider a third electronic state |/) that is coupled 
non-radiatively to the photo-excited state |e) as in Fig. 1(b), the wave-packet created in the excited state by the 
pump can cross over to |/), leaving a vibrationally coherent product jist . Fig. 1(b) suggests the importance of 
taking a multidimensional view of the problem, whereby the surface crossing between the reactant excited state 
|e) and the product state, |/) along the reaction coordinate R is accompanied by the creation of a vibrational 
coherence along the Q coordinate that is coupled to the non-radiative transition. In earlier work, we have 
presented expressions for the time dependent population and first moment evolution of vibrational dynamics 
following a Landau-Zener surface crossing ||2l| ]. These expressions are rigorously valid, provided the quantum 
yield for the reaction is unity (a condition that holds for NO and CO photolysis from heme proteins p^-|2^). 

A common theoretical formulation for pump-probe spectroscopy is based on the third order susceptibility x^^^ 
formalism, which provides a unified view of four wave mixing spectroscopies p5|-p8| with different combinations 
of fields, irrespective of whether they are continuous wave or pulsed. However, the separation of the pump and 
probe events is not clear in this formalism, since the pump induced density matrix is implicitly contained in 
the third order response functions. Thus, it is also attractive to treat the pump and probe processes separately 
in the well separated pulse (WSP) limit. For example, the "doorway/ window" picture has been developed |28|| 
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which can be used to represent the pump and probe events as Wigner phase space wave packets. This readily 
enables a semi-classical interpretation of pump-probe experiments [P9|-pl| . Another view of the WSP limit is 
based on the effective linear response approach. In this approach, the pump induced medium is modeled using a 
time-dependent linear susceptibility p^-^,^j27| . This has the appealing aspect that a pump-probe experiment 
is viewed as the non-stationary extension of steady state absorption spectroscopy. 



(a) 



(b) 
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FIG. 1. (a) Schematic of a two electronic level system, with ground (\g)) and excited (|e)) electronic states that are 
coupled radiatively. (b) Multidimensional view of vibrational coherence created by a Landau-Zener surface crossing. 
Surface crossing along the R-coordinate between the electronic states |e) and |/) leads to vibrational coherence along 
the Q coordinate in the product electronic state. The states \g) and |/) are assumed not be radiatively coupled. 



While the above theoretical treatments are primarily concerned with field driven processes, the need to in- 
corporate additional non-radiative mechanisms into the theory is evident from recent experimental observations 
| p^ , pT|j35| . Recent studies have addressed population transfer during non-adiabatic transitions in electron trans- 
fer systems . Numerical wave packet propagation techniques have also been presented to model coherent 
dynamics due to non-adiabatic excited state processes p8| , p9[ . 

The present work is motivated by the need to interpret FCS experiments of reactive and non-reactive samples 
of myoglobin |3^. Here, we make the WSP approximation, and demonstrate the usefulness and accuracy of the 
effective linear response approach, with the non-stationary states represented using unitary displacement and 
squeezing operators p9| . This representation readily leads to rigorous analytic expressions for the effective linear 
response functions. The analytic expressions are physically intuitive and allow for the efficient calculation of 
oscillatory amplitude and phase profiles over a wide range of carrier frequencies. The resulting information is the 
FCS analogue of a resonance Raman excitation profile. The computational ease of this approach eliminates the 
need for impulsive or semi-classical approximations ||39|-|45|. Furthermore, when the linear response functions are 
used in combination with a moment analysis of the pump induced (doorway) state pq] (presented in Appendix 
A), excellent agreement is found with the full third order response approach l27|j28|| . In addition, coherence 



2 



driven by rapid non-radiative processes as in Fig. 1(b) can be readily incorporated into this formalism. This is 
illustrated using the first moments (presented in Appendix B) of a Landau-Zener driven coherence 

The general outline of the paper is as follows. In Sec. II, we briefly review the third order response approach 
and its connection with the effective linear response approach in the well separated pulse limit. In Sec. Ill, we 
obtain analytic expressions for the effective linear response function using a displaced thermal state representa- 
tion for the doorway state. We further consider undamped vibrational motion and derive arbitrary temperature 
expressions for the dispersed and open band pump-probe signals. In Sec. IV, we present simulations to demon- 
strate the accuracy and physically intuitive aspects of the first moment based effective linear response approach. 
Next, we include the effects of inhomogeneous broadening of the electronic lineshape. Finally, we consider the 
application of the effective linear response approach to reaction driven coherence, and discuss the magnitude of 
overtone signals using displaced and squeezed initial states. 



II. THIRD ORDER AND EFFECTIVE LINEAR RESPONSE 



Here, we briefly review the connection between the third order and effective linear response functions in the 
well separated pulse limit. In the electric dipole approximation, we ignore the spatial dependence of the electric 
fields and write the total electric field of the incident pulses as E(t) — Ea(i) + Eb(i — r). The subscripts 'a' and 
'6' refer to the pump and the probe, and r is the delay between the two pulses. The detected signal in the open 
band detection scheme, is the differential probe energy transmitted through the medium 

r + co 

Sir) = - / dtEb{t - r) • idP{t, T)/dt), (1) 



where F{t, t) is the material polarization induced by the pump and probe fields. If we assume that the delay r is 
larger than the pulse durations, the lowest order in the polarization that is detected in pump-probe spectroscopy 
is given by 

m,T)^ f dh j dh f dhx%{tMMM)Ebj{h-T)Eak{t2)Eai{ti). (2) 

J ~oo J —oo J —oo 

Here, Xi%i is the tensorial third order susceptibility. We have written the sequential contribution ||2^] ignoring 
the coherent coupling/ tunneling terms [|l^,^ that involve the overlap of the pump and probe fields. The 
state of the system in between the pump and probe pulses is not exposed in the the above expression for the 
polarization. If we group the ti and t2 integrations in Eq. (|^) that involve only the pump interactions, then we 
are left with the ^3 integral over the probe field and we may rewrite Eq. (||) as 

m,T)^ [ dt3X-^fit,t3)Etj{t3~T), (3) 
^ — oc 

where xlj^ {^^ ^3) is the effective linear susceptibility describing the non-stationary medium created by the pump 

interaction |2^,^. Recall for a stationary medium that xtj'^ {t, ^3) = Xij(^~^3)j which describes the equilibrium 
response of the system. 

Eqs. (|^) and (^ are formally equivalent expressions for the polarization induced by the well separated 
pump and probe pulses. However, Eq. (^) is more important than a mere rephrasing of the third order 
polarization expression in the well separated limit. The non-stationary response function Xi/'^(^: ^a) can equally 
well describe coherence induced by mechanisms other than (and including) the pump field interaction; e.g. a 
rapid non-radiative surface crossing. To see what is involved, we evaluate the induced polarization quantum 
mechanically (in the interaction picture [|^) as 

P(t)-Tr(p(i,r)/2(i)), (4) 

p{t, t) is density operator describing the quantum-statistical state of the system, and p, is the electric dipole 
moment operator. We write the total Hamiltonian for the molecule and its interaction with the laser fields as 
H{t) = Ho + Hi{t), where the free Hamiltonian Ho and the interaction Hamiltonian Hj have the following form 
for a molecular system with two electronic levels: 
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where Hg and = -ffg + V , are, respectively, the Born-Oppenheimer Hamiltonians for the ground and excited 
electronic states. V is defined as the difference potential that specifies the electron nuclear coupling and Tiil^ is 
the vertical electronic energy gap at the equilibrium position of the ground state. With the above definitions, 
the third order susceptibility takes the standard form 

xfliitMMM) = {i/hfTr{%{t)[ri,{h),[rLk{t2)Ah{ti).K-^)]]])- (6) 

Here, p{—oo) is the initial density operator of the system before the pump and the probe pulses and Tr denotes 
the trace operation. Expansion of the above triple commutator leads to four non-linear response functions and 
their complex conjugates (|J2^. On the other hand, if we substitute Eq. in Eq. (Q), and group the pump 
field interactions separately, a comparison of the resulting expression with Eq. (|^) leads to 

Xl^^{t,h) - {t/n)Tr{Mt) UlAt3),Sp{t3)]) = {i/n)Tr{5p{h) [Ut).pAh)])- (7) 

Here, 5p{t) is the perturbation in the density matrix due to the pump interaction (also referred to as the density 
matrix jump p8[ | and the doorway function p9| , pl[ |) and is given (in the interaction picture) by, 

^P(*3)=('r) f'dt2f dti[Mt2),[Mtl),pi~^)]]Eakit2)Ealih). (8) 
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In arriving at the second equality in Eq. (^, we have used the invariance of the trace under permutation. In 
the WSP limit, where the overlap between the pump and probe pulse is negligible, we can extend the upper 
limit ^3 in Eq. (^ to infinity. The interaction picture density matrix is then time independent. 

Eq. (|^) expresses the effective linear response as the average value of the standard linear response Green's 
function p9[ with respect to the pump induced (non-stationary) density matrix Sp. Note that Sp does not 
necessarily have to be induced by the pump fields as in Eq. (^). It could, for example, describe the perturbation 
in the state of the system due to a non-radiative interaction that follows the pump excitation to a dissociative 
state. It might even be a pump induced state with higher order pump interactions included as would be required 
for very strong pump fields [ pO| . The only assumptions are that the probe field interrogates the non-stationary 
medium after the coherence has been created, and the probe interaction with the medium is linear. 

We express the full perturbed density matrix Sp in terms of the ground and excited electronic state nuclear 
sub-matrices. The point is that the second order density matrix in Eq. (|^) has no electronic coherence due to 
the even number of dipole interactions. Hence 

Sp = Sp,\e)(e\+Spg\g){g\^(^^f . (9) 

Correspondingly, the effective linear susceptibility Eq. (j^) can be decomposed into a ground and excited state 
linear response function 

X^'(t,i3)=x!f(t,t3)+xl;^(i,t3), (10) 

where 

xg'^'(t,i3) = m)Tr [Spe.^g {e,g\ [^.(i), ^.(ts)] \e,g)] . (11) 

In general, the pulse induced nuclear density matrices Spg and Spe contain vibrational coherence, i.e. off- 
diagonal elements in the phonon number state representation. This coherence translates into time dependent 
wave-packets in a semi-classical phase space {Q,P) Wigner representation of the density matrix | |3l| , pl| , p2| , p8[ . 
On the other hand, the highly localized (in Q and P) nature of the impulsively excited non-stationary states 
suggests that we calculate their moments using 

X,,g{t) = Tr [(e, g\ Xit) |e, g) 5% J . (12) 



4 



With a^a) as the creation (destruction) operator for the phonon mode, X represents the dimensionless quadra- 
ture operators Q = (a^ + a)/V2, P = z(a^' — a)/-\/2 or their higher powers. The time t is larger than the pulse 
duration. Let the initial {t = 0) values of the first moments of the position Qg{0) and momentum PsiO), for 
(say) s — e or s — g he denoted by (QsOi Pso)- These uniquely determine the subsequent first moment dynamics 
Qs{t). Since Qso and Pso denote shifts from thermal equilibrium, we may represent the pump induced nuclear 
density matrix for the electronic state s 

Sp,=D{\s)p^^^D\Xs), (13) 
where is the equilibrium thermal density matrix corresponding to the electronic level s. 

p^r^^ = Z-^ exp (-HjkBT^ . (14) 

D{\s) is the quantum mechanical displacement operator defined as 

D{\s) = exp {Ka^ - A*S) . (15) 

Here = (QsO + iPso) /a/S is the initial displacement (in phase space) of the coherent state induced on the 
potential surface s. When Eq. ( p^ ) is substituted into Eq. (pi]), the response functions x''^'^'* are readily 
calculated as we show below. Furthermore, Eq. ( [l^ offers a general scheme to represent the non-stationary 
density matrix on any given electronic level with only a knowledge of the first moment A^. The effective 
linear response function is thus not restricted to coherences driven by a second order pump interaction. This 



is the point of departure of the present work from earlier treatments 33,^,|27|| which generally utilize the 
second order pump induced density matrix Eq. (H), and are thus formally identical to the x^^^ approach. 
This development can be extended to include the higher moments of Q and P, but for harmonic potentials we 
expect the higher moments to play a less significant role in the overall dynamics. Second moment changes are 
incorporated in a manner analogous to Eq. (031) in Appendix D. 



A. Detection schemes 



The two common experimental detection schemes are the open band and the dispersed probe configurations. 
For simplicity, we assume that the medium is isotropic and take the pulse fields to be scalar quantities. We 
then write the susceptibilities without tensor subscripts. The measured open band signal is given by Eq. (|^). 
For an almost monochromatic laser pulse, we can write the electric field and the induced polarization as: 



■fc(f-r) = 5Re|£(i-r)e-*'^<=(*-^)} and F(i, r) = SRe |p(t, r)e-'"<=(*-^)} , (16) 



where ujc is the carrier frequency, and £{t — r) and 'P{t, r) are slowly varying envelope functions. Typically, the 
envelope function for the probe field is a real Gaussian centered at time r; £{t — r) = EQG{t — r), where Eq is 
the field strength of the probe pulse. The envelope function for the polarization can be obtained by using the 
definitions in Eq. ( p6| ) and Eq. (^. After making a change of variables s = i — ^3, we find 

/•OO 

r{t,T)= dse^^^'x^^Ht^t- s)£{t- s-t). (17) 
Jo 

Employing Eq. (|l^), Eq. (Q), and the rotating wave approximation (RWA) in which highly oscillating non- 
resonant terms are ignored, we get 

/oc 
dt£*{t-T)V{t,T). (18) 
-00 

Since for the present case £{t) is real, the measured (dichroic) signal is directly related to 3m{'P(t, r)}. If 
£{t) were made imaginary as in heterodyne techniques [ p4| , pl| , then the resulting (birefringent) signal is related 
^e{V{t,T)}. 

In the dispersed probe detection scheme, the measured quantity is spectral density 5(w,t) of the transmitted 
energy through the sample defined by 
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POO 

S{t) = I S{uj,T)duj. 
Jo 

In the RWA, 

S{uj,t) = -(w/47r)3m {e-^('^-'^<=)^£*(u; - ujc)r{uj - oj^, t)} , 
where f (w) and 1^(1^) denote the Fourier transforms of the envefope functions, £{t) and V{t,T). 



(19) 



(20) 



III. EVALUATION OF THE RESPONSE FUNCTIONS 

In this section, we evaluate both the third order and the effective hnear response functions. Consider the 
ground state response function Xg(^: ^3) given by Eq. (|ll|) without the tensor subscripts. We make the Condon 
approximation and ignore the coordinate dependence of the dipole moment operator. If we then let fl = 
fige \g) (e| + fieg |e) {g\ , and use the interaction picture time evolution /i(i) = e*^''*/''/2e~*^°*''^, we can write, 



where we have defined the two-time correlation function for the ground state response: 



Spg exp 



dsV{s) 



(21) 



(22) 



Here, the subscript + denotes time ordering and V{s) evolves in time s via Hg. The density matrix Spg is 
obtained from Eq. (||) as {g\ Sp\g), assuming p{—oo) — p^^ \g) {g\: 

|2 noG nt2 



6pg 



dt2 



dtiEa{t2)Ea{tl) 



X < e 



-iO„(t2-ti) 



exp -- 



dsV{s) ) p^/^ + h.c 



(23) 



where h.c. denotes Hermitian conjugate. 

For the excited state response function, we similarly write 

where Ce is the two-time correlation function for the excited state response: 



Tr 



(5peexp(^ / (isi/(s)^ exp(^-^ [ dsV{s] 



(24) 



(25) 



/o / - \ " ^0 / + 

with the subscript — denoting anti-time ordering. The excited state density matrix 6pe is obtained as {e\dp\e): 



Spe 



dh 



X < e*^^""* exp I - 



dhEaiti)Ea{t2) 

) 

'-j\s'V{s')^ pi^'>exp(^~^j\s'V{s')] +h.c. 



(26) 



The above expressions are valid for a two level system with arbitrary difference potentials. In what follows, 
we take V — — (?i/ma;o)^^^/Q = — {huJoA) Q, with dimensionless Q and relative displacement A. The electron- 
nuclear coupling force is expressed as f — A(r7ia;Q?i)^/^, where m and luq are, respectively, the reduced mass 
and frequency of the mode. 
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If the trace in Eqs. ( p2| ) and (|2^) were evaluated with respect to the equihbrium thermal density matrices 
yo^'®^ rather than Spe^g, then the correlation functions, Ce,g{t,t^) = K^^git — t^) would be the equilibrium 
absorption and emission correlators [^6| ; i.e. 

Kg{s) = e-^i^oos~r,\,\^-g{s). ^^(^^^ ^ g»Ooos-re|s|g-s(s)^ ^27) 

where we have introduced homogeneous dephasing through the electronic damping constant Fg. For linearly 
displaced and undamped oscillators, g{s) is given by 

g{s) = (AV2) [{2nT + !)(!- cos(t^os)) + i sin(tJos)] , (28) 

where 2nT + 1 = coth (hujo/ 2k bT). Exact expressions for the damped harmonic oscillator |^ can be easily 
incorporated into the present development. The half- Fourier transforms of Kf,^g{t) determine complex lineshape 
f mictions, 

^{uj)=i dse"^'Kg{s), and 6(0)) = i / dse"^' K*{s), (29) 
Jo Jo 

whose imaginary parts are directly related to the absorption and emission cross-section: (7a{lo) cx aj$/(w) and 



A. Full third-order response 

Analytic expressions have been derived for the third order susceptibility x*^'^^ for a two level system coupled 
to a multimode set of linearly displaced harmonic oscillators, and expressed in terms of non-linear response 
functions Rj (j — 1..4) |2^,^^. For completeness, these are listed in Appendix C. When the second order pump 
induced density matrices in Eqs. (p3|) and (p6|) are substituted in Eqs. (E^) and (p5|), the ground and excited 



s tat e correlation functions Cg.e a-re related to the non-linear response functions as expressed in Eqs. (CI) and 



(C4) 



The expressions of Appendix C cast the effective linear response functions in terms of the conventional non- 
linear response functions. The two approaches are entirely equivalent at this stage. The key point from Eqs. 



(CI) and (C4) is that the ground and excited state correlation functions Cg,e{t,t3) involve a double integration 
over the pump electric fields. In the next section, we will show that the displaced state representation for the 
non-stationary density matrix Eq. (Hq) directly leads to analytic expressions for Cg^eitjt^). 



B. Effective linear response 

1. Displaced thermal state 
Following Eq. (|l3|) , we represent the pump induced non-stationary state as 
Sp^D{Xg)pf^'>D\Xg) \g) {g\ + b{K)i^f^ D\X,) |e) (e| 



(30) 



This representation can be used to calculate the effective linear response in a straightforward manner |58| . 
Substituting the ground state part of Eq. (^0|) into the correlation function in Eq. ( p^ and using the permutation 
invariance of the trace and the unitarity of the displacement operator, we obtain 



p^/^ exp ( itJoA / dsD\\g)Q{s)D^ {\g 



(31) 



The effect of the displacement operator is to shift the position operator by a time dependent classical function 
(c-number), namely the mean pump induced displacement 



D\\g)Q{s)D{\g)=Q{s)+Q(s), 



(32) 
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where Q„{s) = ^e{\ge ^"o'*} is the time dependent mean position of the oscillator in the ground state. If 



we substitute the above expression into Eq. (31), then the c-number Q„{s) can be removed outside the time 



ordering and the thermal average simply reduces to the equilibrium absorption correlator Kg(t — t^) in Eq. (|27|). 



Using Eq. (A3) of Appendix A for the first moment dynamics, we find 



Cg{tM) = Kg{t-t:,)e^^ 



iuj^AglS. 



dse cos(wos + ^fig] 



(33) 



where Ag = |Ag| and ipg = arg(A*) = — tan~^ {Pgo/Qgo) are the amplitude and phase for the coherent wave- 
packet motion in the ground electronic state. The above expression has an intuitively appealing form: the 
ground state correlation function for the pump induced non-stationary medium is expressed as a modulation 
of the standard equilibrium (linear) absorption correlator Kg{t ~ ^s) by the first moment dynamics of 
the ground state wave-packet motion. It is clear that the corresponding non-stationary ground state response 
function %g(i,t3), Eq. (pT|), would translate in the frequency domain into dynamic absorption and dispersion 
lineshapes p^,p|,p7| , which determine the final probe response to the non-stationary medium. After performing 
the integral, Eq. (|33|) reads 



Cg(i, ts) = ifg(t - is) exp lAgA/uo^e "''^ {LOySm^cUyt + ipg) ~ jcosiLu^t + ipg)) 
-e^'^*^ (cj„ sin(a;„t3 + ipg) - j cos{uj.vt3 + ipg)) | . 



(34) 



It is interesting to note that the strength of the first moment modulation appears through the product of the 
initial displacement and the optical coupling, Ag/S.. 

For the excited state correlation function, we similarly find 



Ce(t, is) = ■f«^e(t-i3)exp - lAe/^/ujo^e {ujyiiin{ujyt) - ^cos{ujyt)) 
-e"''*^ (w^ sin(cjt,i3) - 7COs(a;„t3)) | . 



(35) 



Thus, the excited state non-equilibrium correlation function is expressed as the modulation of the equilibrium 

(fluorescence) correlation function by the excited state wave-packet. 

The analytic expressions in Eqs. (^) and ( ^5| ) together with the first moments presented in Ens. (Al), 
( [A^ ) and (A4), effectively replace the non- linear response expressions derived in Eqs. (CI) and (C4). The 



generalization of the above single mode results to the multimode case is straightforward. It easy to show that 
the multimode expressions for Cg and Ce, factor into a product of single mode correlation functions. The results 
derived using the first moments are approximate, since the higher moments of the wave-packet motion (which 
are induced by the pump pulse [4^]) are neglected. Nevertheless, they provide a fully quantum mechanical 
description of the probe response due to modulation by coherent nuclear dynamics. Moreover, we expect that 
for a harmonic system, the first moment dynamics will constitute a major part of the wave-packet motion 
detected by the probe. This expectation is verified when we present simulations comparing the first moment 
based linear response approach with the full third order approach. 



2. Analytic expressions for pump-probe signal: CARS and CSRS responses 

An advantage of the analytic expressions in Eqs. (|^ and (^5|), is that for the special case of undamped 
nuclear motion, they can be expanded in a Fourier-Bessel series, from which the polarization Eq. (p7| ) can 
be evaluated in closed form. This yields an analytic expression for the dispersed probe signal Eq. (eM that 
provides insight into the origin of both resonant and non-resonant pump-probe signals. If the series expansion 

oo 

e'^^'"^'^^ J2 exp{im9)Jm{g), (36) 

m— — oo 

is placed into Eq. (jsj), we find for 7 = that 

CSO CO 

Cg{t,t3)^Kg{t~t3) ^ e^M-o*3+^.) ^ J„(AgA)J„_fc(A<,A)e™"°(*-*^). (37) 

k—~oo n— — oo 
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We put this expression into Eqs. (^T|) and to find the polarization within the RWA (which neglects the 
contribution of C*). With Gaussian pulses, S{ijj) = EoG{uj), where G{u!) is a Gaussian spectral function. The 
ground state contribution to the frequency dispersed pump-probe signal is finally obtained: 



Sg{uj,T) = ^Cfe(a;)cos(fc (wot + (pg)) + Sk{uj)s'm{k {ujqt + ipg)) 

k=0 

oo 

fc=0 

where the amplitude and phase of the fc'th overtone are 

= kipg - tan"^ [Sk{uj)/Ck{uj)] , 
and where Ck{uj) and Sk{uj) are the quadrature amplitudes given by, 

SkiL^)^coM\Gp{LOd,kuo)<p'-^\Lo)~GpiLOd,~ku;o)<i>^-''\uj) 



(38a) 
(38b) 



(39a) 
(39b) 



(40a) 
(40b) 



In the above expressions, we have defined the constant M = 
the product spectral function of the probe pulse: 



'{\fige\'^EQ/4:TTfi) for convenience, and introduced 



Gp{ujd, kujo) = G{uJd)G{uJd + kuo), (41) 
where ujd ~ lo — ujc the detuning frequency and the absorptive and dispersive basis functions are defined as 



$f'(w)= JniAgA)Jn^kiAgA)^iicj + nu;o), 



n — — OQ 



$^*^'(w)= J2 JniAgA)Jn^kiAgA)^Riu; + nu;o). 



(42a) 
(42b) 



The basis functions depend only on Ag, A and the equilibrium lineshape functions and do not depend on 
the properties of the probe pulse. Note the symmetry conditions $j ''\<-^) — ^^'^^lo — kuJo) and $^ '"''(f^) = 

^'r\uj - kujo). Eqs. (38-42) are rigorous expressions for the dispersed probe signal for a single undamped 
mode coupled to a two electronic level system. They are valid for arbitrary temperature, detection frequency 
and pulse width. The simple and direct connection between pump-probe signals and the equilibrium lineshape 
functions is exposed by these results. Equilibrium lineshapes thus enter the calculation at t wo separa te s tages: 
firstly in the expressions for the pump induced first moment amplitude and phase in Eqs. ( |AlD and (A2), and 
secondly in the final probe detected signal in Eqs. (42a - 42b ). 

It is clear from Eqs. (38) that the detuning dependence of the dispersed probe signal is determined by the 
amplitudes Cfc(w) and Skit-^), which are in turn related to the product of the spectral function Gp and the 

displaced lineshape functions ^fli- In studying the detuning dependence of Sg{uj,T), it is interesting to look 
at two opposite limits. When the electronic dephasing time is much longer than the pulse durations and the 
vibrational periods, the basis functions $ 
that acts like a filter in Eqs. (40a) and 
the electronic dephasing time is much shorter than the pulse durations and vibrational periods. In this case, 
the lineshape functions are broad, exhibiting a much slower variation with respect to u; than the pulse envelope 
spectral function which then acts like a filter in Eqs. (40a) and ([40b| ). The dispersed signal then consists of 
a superposition of red and blue shifted field envelope functions Gp(wd, ifci^o)- The superposition weighting 
depends on the amplitude of the initial pump induced displacement and the displaced equilibrium lineshape 
functions. Thus, resonances occur at frequencies to ~ uJcTkoJo/'^ which correspond to the peaks of Gpiuid, ifcwo). 



'lo) and (w) consist of a well resolved Franck-Condon progression 
Oq) and determines the detuning dependence. In the opposite limit. 
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These resonances can be identified with the weh known coherent Stokes Raman scattering (CSRS) and coherent 
anti-Stokes Raman scattering (CARS) signals |60|j6l|] . 

The delay dependence of the dispersed signal is composed of oscillations at all harmonics of the fundamental 
frequency wq- The amplitude and the phase of the /c*'* harmonic is given by Ak{u}) and (/>fe{a;). It is clear 
from Eq. ( ^9b| ) that the dispersed probe signal phase 0fe(w) is not simply related to the initial phase ipg of the 
coherent motion because of the additional frequency dependent functions Ckiuj) and Sk{to). However, it can be 
easily shown that when we integrate Sg{Lu,T) as specified in Eq. (p^, the integral over Sk{uj) is sm all d ue to 
approximately cancelling contributions from the Stokes and anti-Stokes shifted components in Eq. ( [40b| ) . We 
then get for the ground state contribution to the open band signal 

oo 

Sgir) =^ Ask COS {kujoT + (l)sk) , (43) 

k=0 

where the amplitude and phase of the fc'th overtone are given by 

/•oo 

Ask - / dLj\Ck{u;)\, (44a) 
Jo 

0sfc = fcf^g - [0 or tt]. (44b) 



Thus, the phase of the fundamental open band signal (f>si determines the initial phase cpg of the nuclear motion 
to within an additive constant of tt. The additive constant phase arises from the fact that Cfc(tj) can have 
a positive or negative value depending on whether the initial transmission of the probe pulse is increased or 
decreased. It is easily shown that the open band signal will be of the form in Eq. (^3|) even when several 
vibrational modes are active. These results demonstrate that the phases of the open band signals are a direct 
reflection of the initial conditions of the non-stationary states. 

It is also notable from Eqs. (43-44) that the open band (dichroic) signal vanishes as we tune off-resonance p2| , 
where the imaginary lineshape function is negligible and Ck{uj) is vanishingly small. In contrast, the off-resonant 
frequency dispersed signal is non- vanishing and depends on the real part of the lineshape through Sk{Lo). We 
also have (pg —11/2 off-resonance using the results of Appendix A. We then get 

00 

§<,(o^,t) =^5fc(u;)sin(fc(t^oT + 7r/2)). (45) 

k=0 

The dispersive functions ^'"^\uj) are approximately constant in the off-resonant limit. It then follows from Eq. 



(401:) that the detuning dependence is mainly determined by the pulse spectral function Gp. The dispersed 
probe signal consists of CARS and CSRS resonances that are centered near u ^ uicT and oppositely 

phased 

The expressions derived here make minimal assumptions and can be readily extended to incorporate optical 
heterodyne detection schemes p] , ^2|j5^ , |55| ] . We have derived explicit expressions for the dichroic response 
Sm{'P(i,T)}. It ca n b e shown for the birefringent response 5Re{7'(i, r)}, that the roles of the quadrature 



amplitudes in Eq. (38a) will be reversed; i.e. that Ck{Lo) is the coefficient of the sine term and Sk{Lo) is the 
coefficient of the cosine term. Also, the CARS and CSRS contributions of Sk{uj) will add constructively to 
give a non-vanishing birefringent open band signal off-resonance. This has been observed in transparent liquids 
p5[ . For the resonant case, the results derived here are a multilevel, high temperature generalization of earlier 
treatments based on density matrix pathways which considered a pair of vibrational levels in the ground 
and excited states. In contrast to prior treatments no assumption has been made with regard to the 

pulse durations in deriving Eqs. (38-42). Finally, we note that although we have given expressions for only the 
ground state response, analogous results are easily obtained for the excited state response. Here, the emission 
lineshape function Q{uj) in Eq. (E9|) plays a role similar to $(a;) in the ground state expressions. 



3. Reaction driven coherence 



As previously discussed, an advantage of the effective linear response approach is that it allows for a rigorous 
calculation of the probe response to non-radiatively driven coherence. As an example, we consider the multilevel 
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system depicted in Fig. 1(b), and apply the effective linear response approach to the detection of reaction induced 
coherent nuclear motion along the Q degree of freedom in the product state |/). 

The chemical reaction step and subsequent probe interaction minimally constitute a three electronic level 
problem comprising the electronic states |e), |/) and |/'). The states |e) and |/) arc non-radiatively coupled. 
The probe interaction couples the ground state |/) and excited state |/') of the product, and is assumed to be 
well separated from the reaction step. The Hamiltonian for the reaction-probe stage of the problem may be 
written as 

/ Hf, -Hf'Eb{t) 
HNRit)^ \-^if,fEt,{t) Hf + UfiR) J |. (46) 

\ J 

Here, Uej are the dissociative potentials along the classical reaction coordinate R, for the states |e) and |/). Hej 
represent the quantum degree of freedom (Q) coupled to the non-radiative transition. J is the non-radiative 
coupling parameter which, for the case of MbNO photolysis discussed below, corresponds to the spin-orbit 
coupling operator needed to account for the spin change of the heme iron atom upon ligand photolysis. For 
simplicity, we also assume that the pump interaction that couples the states \g) and |e) is well separated from 
the chemical reaction. The sequence of interactions and the relevant coupling constants is then: pump, fiEa 
chemical reaction, J —^ and probe, fiEi,. Each individual interaction is assumed to occur independently. The 
separation of these events allows us to consider the four electronic level problem in Fig. 1 (b) as three sequential 
two electronic level problems. 

We make the following two assumptions which are satisfied in the case of MbNO photolysis: (i) Before the 
curve crossing takes place, the system (in reactant state |e)) is assumed be in thermal equilibrium along the 
Q coordinate, (ii) The quantum yield for the dissociative reaction along R is assumed to be unity; i.e. all the 
reactant molecules are transferred to the product state during the surface crossing. Assumption (i) holds if the 
Q coordinate is not optically coupled to the \g) and |e) electronic states. The pump pulse merely transports 
a fraction of the ground state electronic population to the excited state, leaving the vibrational state along 
Q unchanged. However, the R degree of freedom is dissociative on the excited state potential surface and is 
left far from equilibrium after the pump excitation. For example, the 220 cm~^ Fe-His mode in MbNO is not 
optically coupled as revealed by its absence in the resonance Raman spectrum of MbNO. Its strong presence in 
the pump-probe signals demonstrates that this mode is triggered into oscillation following the highly efficient 
process |22 -24| of ligand photo-dissociation in MbNO 



A multidimensional Landau-Zener theory was previously developed to describe vibrational coherence induced 
by non-radiative electronic surface crossing based on assumptions (i) and (ii) above. For completeness of 
presentation, we have reproduced the key results of that work in Appendix B. We also present expressions for 
the amplitude Af and phase ff of the first moment of the nuclear motion induced on the product surface 
I/). Using these reaction driven initial conditions, we make the following representation for the non-stationary 
nuclear density matrix in the product state /Oy-; 

p'y = 5(A/)p</'St(A;), (47) 

^( f) . 

where Xf — Afe^^^ is the complex displacement of the reaction induced coherence and pjl is the equilibrium 
thermal density matrix for the nuclear Hamiltonian Hf. Since only the ground state of the product |/) is 
initially populated, the probe response due to the nuclear dynamics in the product state potential well can 
now be obtained along the same lines that led to Eq. (^J). The correlation function C/ for the non-stationary 
response is obtained by replacing Ag, ipg and A by A/, ipf and A/ in Eq. (^). Here A/ is the dimensionless 
coupling associated with the coupling of Q to the f ^ f transition. The equilibrium optical absorption 
correlator is replaced by Kf{t — t^), which refers to the pair of electronic states |/) and |/'). 

IV. SIMULATIONS AND DISCUSSION 

The numerical simulation of the open band signals involves multiple integrals of the time correlation functions 
which can be evaluated using standard algorithms. The full third order response approach involves quadr uple 



integrations over time. The effective linear response approach replaces the correlation functions in Eq. (CI) 
and (C4) by the analytic expressions in Eq. (|^ and (^). For field driven coherences, the analytic expressions 
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require the evaluation of the first moments in Eqs. ( [Ai1 - |A^ ) and (A4) of Appendix A. The moments can be 
evaluated in a single step before doing the time integrals. Hence, two of the time integrations are eliminated 
and the overall computation time is significantly reduced. The signal amplitude and phase in either of these 
approaches are obtained by performing a Fourier transform of the delay dependent pump-probe signal in Eq. 
(p^). An alternative method that can be used when vibrational damping is neglected, is to directly calculate 
the amplitude and phase of the dispersed and open band signals using the analytic expressions derived in Sec . 
III.B.2. For calculations involving only a few modes, the analytic formula for the dispersed signal in Eq. (38b) 
can be easily evaluated, and a subsequent integration over frequencies as in Eq. ([19|) leads to the open band 
amplitude and phase. 

In what follows, we demonstrate that amplitude and phase excitation profiles of field induced coherence are 
equally well predicted by the third order and the effective linear response approaches. We then apply the 
effective linear response approach to calculate coherent signals induced by rapid non-radiative reactions. 



A. Temperature and carrier frequency dependence 



One of the key parameters in a pump-probe experiment is the carrier frequency of the pump and probe 
laser pulses. Not only is this parameter the most easily accessible to the experimentalist, but measurements of 
the oscillatory amplitude and phase profiles through the resonant region provide crucial information regarding 
the origin of vibrational coherence The computational simplicity of the effective linear response 

approach, using the first moments, enables a direct calculation of phase and amplitude profiles over the entire 
absorption spectrum using realistic pulse widths. 

In order to illustrate the accuracy and the physically intuitive aspects of the effective linear response approach, 
we treat a simple model system. Consider a single undamped mode with luq = 40 cm~^ and S = (A^/2) = 0.5 
coupled to a homogeneously broadened two level system with Fe — 800 cm~^. We assume 10 fs pump and probe 
pulses in a degenerate (same color for pump and probe) configuration. The short pulse width allows comparison 
with the third order response calculation which is quite formidable for long pulse durations ~ (50 — 100) fs. 

Expressing the delay dependent oscillatory signal as 

S{t) ^ A, cos(woT -I- <j)s), (48) 

we plot in Fig. 2 the degenerate open band amplitude As and the phase 4>s profiles for the 40 cm~^ oscillations. 
The ground and excited state profiles are plotted for a range of carrier frequencies across the absorption maxi- 
mum. The top panels of Fig. 2(a) and (b) show the ground and excited state amplitude and phase profiles for 
T ^ K and T = 300 K. Both the third order response and the effective linear response outputs are plotted. 
The excellent agreement between the two approaches is evident over the entire range of carrier frequencies. 
We note that the ground and excited state amplitudes dip near the classical absorption and emission peaks, 
respectively, fly and fl^ — 2liJqS. The phase of the ground state signal shows a variation of 27r as the carrier 
frequency is detuned across fi^. There is a sharp jump between ±(7r/2) at r2„. In contrast to the ground state, 
the excited state phase remains constant apart from a tt phase jump at ri^, — 2aJo*S'- At low temperature, the 
amplitude of the ground state signal drops by almost an order of magnitude. The approach of the phase toward 
±(7r/2) on either side of the discontinuity at is steady and almost linear. For the high temperature case, 
the approach of the phase toward ±(7r/2) occurs much more sharply near Vly. While the excited state phase 
is independent of temperature, the excited state amplitude becomes more asymmetric as the temperature is 
increased, as seen in Fig. 2(b). 

All of the above aspects of the ground and excited state signals can be clearly understood using the first 
moments of the pump induced oscillations; see Appendix A. First, recall that the open band phase yields the 



initial phase of the pump induced oscillation to within an additive factor of tt as in Eq. (441 ). This is shown for 
the ground state signal in Fig. 3. We plot the open band phase of the 40 cm~^ mode over an expanded range 



that includes the off-resonant limit, along with the initial phase of the wave-packet calculated using Eqs. (jAlJ) 



and (A2). 
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T= K Ground state T= 300 K 





FIG. 2. Comparison of one-color pump-probe open band amplitude As and phase <j>a profiles predicted by the first 
moment based effective linear response approach (solid line) and the third order response approach (circles) . The profiles 
shown are for a single undamped mode {luo = 40 an~^ , S — 0.5) coupled to a homogeneously broadened (Fe = 800 cm~^) 
two level system at T = 300 K and T — K. The top panels of (a) and (b) show the amplitude and phase profiles for 
pulse width of 10 fs. The bottom panels of (a) and (b) show the amplitude of the ground and excited state signals as 
two dimensional functions of the pulse width Tp and the pulse carrier frequency ujc, for both T = K and T = 300 K. 
The analytic expressions in Eq. (38-42) were employed in the 2-D simulations. 
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Fig. 3(c) schematically depicts the effective initial conditions prepared by the pump pulse in the ground and 
excited states. Shown are the initial conditions for three different pump carrier frequencies, near resonance 
(indicated by downward pointing arrows in Fig. 3(a)) and off-resonance towards the red side of the absorption 
maximum. The arro ws abov e the wave-packets indicate the direction of the momentum induced by the pump 
pulse. From Eqs. (A1-A2), we note that the mean nuclear position and momentum in the ground state 
depend on the derivative of the absorption and dispersion lineshapes respectively [Q. Thus, the wave-packet 
is undisplaced from equilibrium for excitation at loc — The momentum attains a maximum value at this 
frequency, and is signed opposite to the excited state equilibrium position shift j4^E3,H. 
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FIG. 3. Interpretation of the open band phase for the ground state signal in terms of the first moments of the pump 
induced nuclear wavepackets. (a) The phase (j)s of the 40 cm~^ mode (from the example in Fig. 2) at T = iC' is 
plotted (solid line) along with the initial phase tpg of the pump induced wavepacket (circles), (b) Same as in (a) but 
for T = 300 K. (c) Schematic of the pump induced efi^ective initial conditions in the ground and excited states in a 
two level system, shown for three different pump/probe carrier frequencies across the absorption maximum n„, and for 
non-resonant excitation ujnr- The arrows above the wavepackets indicate the direction of the momentum imparted. The 
panels directly below the potential curves show the corresponding oscillatory signal as a function of pump-probe delay 
time. 



As depicted in Fig. 3, when the pump pulse carrier frequency is tuned toward the red side of the absorp- 
tion maximum w^, the ground state wave-packet is displaced toward decreased bond-lengths and receives a 
momentum kick in the same direction. Thus, a probe pulse of carrier frequency lor incident on the sample 
sees an initially depleted nuclear distribution (bleach) in the ground state. The probe difference transmission 
signal (pump-on minus pump-off) will at r = be positive and will simultaneously have a positive slope due 
to the sign of the wave-packet momentum. This is seen from the delay dependent signal plotted directly below 
the potentials in Fig. 3(c). The phase of the oscillatory signal S{t) in Eq. (|4|) will obey (37r/2) < cj)s < 27r. 
On the other hand, for blue detuning from absorption maximum ujb, the initial transmission signal will again 
be positive due to a depleted nuclear distribution in the probed region. However, the signal has a negative 
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slope at r = 0, since now the momentum of the wave-packet points toward the probing region. In this case 
the signal phase will obey < (j)s < {tt/2) and is also precisely in phase with the wave-packet itself. When the 
carrier frequency is tuned to resonance at tOy, the pump pulse merely imparts a momentum to the wave-packet 
so that the wave-packet phase is precisely (7r/2). As described below, the CSRS and CARS resonances are 
oppositely phased for excitation at resonant maximum. The integrated signal therefore vanishes at lu^ and the 
phase becomes undefined. Incrementally to the red and blue sides of the absorption maximum, the open band 
phase shows a discontinuous jump of tt since the momentum impulse always points in the same direction in 
the resonance region. Note that Fig. 3(c) shows, for ujy, the oscillatory signal incrementally to the red of the 
absorption maximum, before the phase jump. 

A similar analysis is possible for the excited state signal. The excited state response arises from the stimulated 
emission (by the probe) from the vibrationally coherent excited state. Since the excited state wave-packet 
oscillates about the shifted equilibrium position A, the amplitude dip and phase flip of the pump-probe signal 
from the excited state occurs near the classical emission peak at flv — 2uJoS. As depicted in Fig. 3(c), the excited 
state wave-packet does not receive a momentum impulse; see Appendix A. Furthermore, it is created on the 
same side of the excited state harmonic well irrespective of the pump carrier frequency |46|| . Thus, the phase 
of the excited state wave-packet is independent of the pump carrier frequency. The corresponding signal phase 
is fixed at for blue detuning where there is increased stimulated emission at t = 0. The phase is fixed at tt 
for red detuning where the stimulated emission is minimum at r = 0. Note in Fig. 2(b) the nearly order of 
magnitude increase in the amplitude of the excited state signal compared to the ground state. The increase can 
be traced to the fact that for impulsive excitation using weak fields, the excited state wave-packet is situated 
very close to the vertical energy gap 51^. This corresponds to the ground state hole that has a displacement 
much smaller than A ||5^,^,Q. Thus, the amplitude of excited state oscillations about the equilibrium position 
A is much larger than that of ground state oscillations about zero. The excited state coherence will thereby 
be the dominant contribution for systems with long-lived excited states 11 2^. The asymmetry in the 

excited state amplitude profile can be seen in Fig. 3(c) to arise from the fact that the wave-packet amplitude 
in the excited state is larger for blue excitation than for red |^ . 

The temperature dependence of the oscillatory amplitude and phase in Fig. 2 can be understood via the 
relative temperature dependence of the initial position and momentum induced by the pump pulse. For hig h 
temperatures, the mean thermal phonon population n >> 1 and Qgo is enhanced according to Eq. (Al), 
so that the complex displacement — {Qgo + iPgo) /V^ is dominated by the real part. Hence the initial 
phase of the wavepacket stays closer to zero or vr at higher temperatures, except near fly, where Qgo vanishes. 
Correspondingly, the open band phase also stays closer to zero (or 27r) except for the discontinuous transition 
near fly. When the temperature is lowered, Qgo drops sharply since n << 1 and becomes comparable in 
magnitude to PgQ. Thus, at low temperatures, the amplitude of the ground state signal drops dramatically, 
and the phase varies almost linearly on either side of the discontinuity at fly as in Fig. 2(a). The temperature 
dependence of the excited state first moment Qeo is primarily determined by the relative magnitude of <&/ (w) 
and the difference lineshape A<I>/((jj), which appears as the coefficient of fl in Eq. (A4). The difference lineshape 
is responsible for enhancing the asymmetry in the amplitude profile at high temperatures as in Fig. 2(b). 

An important observation to be made from Fig. 3(c) is that for the linearly displaced oscillator model, the 
amplitude and phase profiles of field driven coherence are independent of the sign of the electron nuclear coupling 
A. If A were negative, then the same arguments given above would apply. The amplitude and phase behavior 
would be the sa me as in Fig. 2. This can be seen directly from the fact that the non-linear response functions in 
Eqs. (|C2a - g2b ) and ( CSa - CSb ) depend only on A^ through the function g{s). However, the schematic depicted 
in Fig. 3 clearly indicates that the absolute sign of the potential displacements are not revealed in a pump-probe 
experiment. 

The depiction of pump induced initial conditions in Fig. 3 contradicts prior predictions arising from time 
dependent wave-packet pictures of impulsive stimulated light scattering |^,||,^ . The prior work suggests that 
the ground state wave-packet is always created on the side of the ground state well that is closer to the excited 
potential minimum. This is not correct, as is shown by a careful analysis of the pump pulse interaction p6| , ^ , ^ . 
It is clear from Fig. 3(c) that the centroid of the ground state wave-packet induced by impulsive excitation is 
strongly sensitive to the carrier frequency of the laser pulse. 

From Figs. 2 and 3 we conclude that the moment analysis based approach offers a physically intuitive and 
accurate interpretation of the observed amplitude and phase behavior of open band FCS signals. Another 
significant advantage offered by this approach is that the analytic form of the effective linear response functions 
reduces the computation times. For the comparisons made in Figs. 2 and 3 (top panels), we chose a very 
short 10 fs pulse in order to make the calculations using the third order approach possible with reasonable 
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computation times. However, in many experimental situations, including those reported here, the pulse width 
is often much longer (~ 100 fs). For long pulse widths the two extra integrations involved in the third order 
approach makes the computation quite formidable. In Fig. 4, we compare the computation times for the linear 
response and third order response calculations as a function of pulse width Tp. The advantage offered by the 
effective linear response approach is clear, especially when it is of interest to calculate the amplitude and phase 
profiles at many carrier frequencies. The computation time for calculating the open band amplitude and phase 
using the analytic formulae in Eqs. (38-42) for a single mode and a two mode case is also plotted in the figure. 
The integral of the analytic expressions is over the spectral profiles of the probe pulse so that a larger value 
of Tp implies a narrower range of integration. In contrast to the numerical integrations in time domain, the 
computation time in this case therefore decreases with Tp. As the number of modes involved increases, the 
summation in Eq. ( p8b| ) involves multiple indices and the calculation becomes more complex. However, the 
usefulness of the analytic expression for a few modes is evident from Fig. 4. 
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FIG. 4. Comparison of computation times for the calculation of the one-color ground state pump-probe signal, using 
the third order response approach and the effective linear response approach, as a function of pulse width. Both the 
absolute time to calculate the signal for 40 u>c points and the ratio of the computation times for the two approaches are 
plotted. Also shown is the computation time for calculating the profiles using the analytic expressions for the pump-probe 
signal in Eqs. (38-42). 



As an illustration of the use of the analytic expressions, we calculate the ground and excited state amplitude 
profiles for a wider range of pulse widths (5 — 500) fs in the bottom panels of Fig. 2(a) and (b). The three- 
dimensional plots are useful in capturing the behavior of pump-probe signals over a wide range of pulse widths 
and carrier frequencies. It is seen that the temperature dependence of the ground state signal amplitude is 
rather uniform over the entire manifold of Tp and lOc- The excited state signal shows more interesting behavior, 
with the asymmetry in the amplitude profiles more pronounced at longer pulse widths. Both the ground and 
excited state amplitudes vanish in the limit of very short and very long (compared to the vibrational period) 
pulses, attaining an optimum for some intermediate value of Tp. The vanishing of the ground state signal for 
Tp ^ can be traced to the fact that the ground coherence vanishes in the short pulse limit [p3|p|j27|,^ . It is 



interesting to note from Eq. (A4) that the excited state nuclei are however coherent for Tp ^ 0. It can be shown 



that the corresponding density matrix is simply thermal density matrix p^** placed vertically above Q = on 



the excited state potential well Eq]. Note that the representation in Eq. (O) is essentially exact in this limit. 
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The vanishing of the excited state open band signal for Tp ^ occurs rather due to the fact that the open band 

signal for a spectrally broad probe is essentially the integral over the derivative-like basis functions, Qf^]^{uj), 

which are analogous to plotted in Fig. 5 (see below). 

Finally, we consider a two color pump-probe experiment, which employs a fixed carrier frequency for the 
pump and a variable carrier frequency of the probe. The analysis presented so far applies equally well to two 
color pump-probe experiments. In this case, the pump imparts a fixed initial condition to the wave-packet that 
is determined from the results of Appendix-A. The signal amplitude and phase profiles as a function of the probe 
carrier frequency will be similar to that of the excited state degenerate pump-probe signal. The amplitude will 
exhibit a dip at Qy, and the signal phase will be the initial wave-packet phase, but for a tt jump at f2y. Since 
the pump induced initial conditions are independent of the probe carrier frequency lUc, the computation of two 
color pump-probe profiles is less cumbersome than the degenerate pump-probe profiles. 



B. Dispersed pump-probe measurements and off-resonant response 

When the pump pulse is tuned off-resonance, A^j{uj) — > and Eqs. ( [Ai1 - [A^ ) show that the ground state 
wave-packet merely receives a momentum kick, and is undisplaced from its equilibrium position. Thus, the 
phase of the open band signal approaches ±(7r/2) as seen from Figs. 3(a) and (b). Furthermore, the dependence 
of the wave-packet momentum on A$fl implies that the momentum impulse changes direction as we begin to 
move off-resonance on the red and blue sides of the absorption maximum. Correspondingly, the signal phase 
crosses 27r (zero) into the first (fourth) quadrant as seen in Fig. 3(a) and (b). The approach of the open band 
phase toward the off-resonant limit of ±(7r/2) is more dramatic at lower temperatures, where Pgo is comparable 
in magnitude to Qgo- At high temperatures, the magnitude of Qgo is dramatically enhanced relative that of 



Qgo because of the thermal factor in Eq. (Al). This enhances the contribution of A$/(a;) to the signal phase. 
The off-resonant limit is thus attained at further detuning from absorption maximum at high temperatures. 

While the open band (dichroic) signal approaches zero off-resonance, we showed earlier that the dispersed 
probe signal does not vanish in this limit and is related to ^r{uj) as in Eq. (p5[). In Fig. 5, we demonstrate 
the on- and off-resonance behavior of the dispersed signal, using the same model parameters as in Fig. 2. 
In the top panel, we plot the basis functions for the T = ground state fundamental (fc = 1) oscillations 
^\^\uj) and ^']^\uj) given by Eqs. (42a- i2b| ). The basis functions for fc = 1 are seen to behave roughly as 



the first derivatives of the absorption and dispersion lineshapes. Similarly, the basis functions for the first 
overtone {to) and (w) can be shown to behave like the second derivatives of the lineshapes and so 
on. The Gaussian spectral profile of the laser pulse (r^ = 50 fs) used in the calculation is also shown in the 
upper panel for three different carrier frequencies Wd, uJc2, and cocs- The functions Ci{(jj) and Si{oj) and the 



dispe rsed probe amplitude and phase Ai{uj) and 0i(w) are plotted directly below. As specified in Eqs. (40a) 



and (40b), the quadrature amplitudes Ci{lu) and Si{u!) are obtained by superposing the Stokes and anti-Stokes 
shifted product functions Gp{LUd, llIo)^\^]^{uj) and Gp{uJd, —loo)^\~^\uj), which correspond to CSRS and CARS 
resonances respectively. While the (dichroic) Stokes and anti-Stokcs components are additive for Ci{uj) they are 
subtractive for Si{u;). It is clear from Fig. 5 that at the absorption maximum ujc — tOci, the derivative nature 

of ^f'\uj) gives rise to Stokes and anti-Stokes components of Ci(w) that are oppositely signed on either side of 
the carrier frequency. Thus, the dispersed signal is tt out of phase for red and blue detuning from the carrier 
frequency. The integral of both Ci (w) and Si (w) vanishes, giving rise to a dip in the open band signal amplitude. 
When the carrier frequency is tuned near the shoulder of the absorption spectrum lOc = ^c2 , the Stokes and 
anti-Stokes components of Ci{uj) add constructively to give a single peak. Si{lo) is very small compared to Ci(w) 
since ^'■i^\uj) are small in this region. For off- resonant excitation ujc — Wc3, the dispersive term Si{uj) begins 
to dominate and the signal approaches the limit in Eq. (^) . While the details of the behavior of the dispersed 
probe signal are highly mode specific |^^, the simulations shown in Fig. 5 illustrate the general aspects of the 
dispersed pump-probe signals based on equilibrium lineshape functions. 

It is clear from Fig. 5 that the direct connection with equilibrium lineshapes made in Eqs. (38-42) allows a clear 
picture of the dispersed probe measurements. Furthermore, when dealing with complex multimode systems, 
one can directly use the experimentally measured absorption lineshape and its Kramer-Kronig transform, the 
dispersive lineshape. The need for modeling the detailed mechanisms responsible for line broadening is thus 
avoided. However, additional steps are required when inhomogeneous broadening is significant (see below). 
Generally, this approach is analogous p% to the use of transform methods in the calculation of resonance 
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Raman scattering cross-sections |]68|- [70| , ^ . 
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FIG. 5. Dispersed pump-probe signal for tlie two electronic level system conside red i n Fig. 2 



F.{1), 



uj) and "I> 



(1)/ 



defined in Eqs. 



the absorptive and dispersive basis functions $y 
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([42a|) and ([l2b|), for T = K. The 



The bottom panels show 
and the amplitude Ai{lli) and phase 01 (tj) of the 40 cm~^ 
fundamental of the dispersed signal. The open band phase for each carrier frequency is also shown (filled circle). 



C. Effects of inhomogeneous broadening 

Until now, we have assumed that the molecular system under consideration has a fixed electronic energy gap 
(rioo) and a homogeneously broadened absorption spectrum. The homogeneous broadening mechanism reflects 
the electronic population decay as well as the rapid fluctuations that take place between the system and the 
environment. In addition to the homogeneous contribution, the width of the absorption lineshape can also have 
a quasi-static origin. For example, a distribution of electronic transition frequencies arising from local changes 
in the environment or disorder in the system fFH^ . Several spectroscopic techniques exist |7^-Q , which can be 
used to selectively probe the effects of inhomogeneous broadening. In this section, we consider the effects of 
inhomogeneous broadening on pump-probe signals. Experimental examples of open band phase measurements 
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on myoglobin (Mb) are presented to demonstrate the presence of an asymmetric inhomogeneously broadened 
absorption spectrum for Mb. 

In the presence of inhomogeneous broadening, the homogeneous pump-probe signal must be averaged over the 
inhomogeneous distribution J-'j{floa). In the third order response approach, this averaging is simply a Fourier 
transform of the function J-j a nd c an be carried out independently of the time integrations over the pump 
fields lH^; see Eqs ( |C8a|) and ( C8b| ) of Appendix-C. In t he e ffective line ar response approach, however, the 
averaging over JF/(rioo) cannot be carried out as in Eqs. ( C8a| ) and (C8b). The time integrations over ti and 
t2 involved in the pump interaction are carried out separately to find the first moments, r^oo appears through 
the lineshape functions <I>/(w) and $_r(w) in the first moment expressions, Eqs. ( |A11 - [A^ ) ) . In order to include 
the inhomogeneities, we first calculate the homogeneous pump-probe signal and then numerically average the 
signal over the distribution J^/(f2oo)- When experimentally measured lineshapes are used in place of the time 
correlator expressions, the homogeneous lineshape must be deconvolved from the measured lineshapes before 
being used in the pump and probe steps of the calculation. The inhomogeneous pump-probe signal is obtained 
by a final convolution of the homogeneous pump-probe signal with the inhomogeneous distribution. 

Inhomogeneous broadening affects the ground and excited state responses in a distinct fashion, owing to their 
contrasting phase behavior. Recall that the oscillatory phase of the ground state signal is approximately in 
phase on either side of the resonant maximum fly. Since the inhomogeneous pump-probe signal is essentially 
a superposition of homogeneous signals, the inhomogeneous amplitude will be non- vanishing at fly The phase 
discontinuity will also be smeared out in the presence of inhomogeneity. 
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FIG. 6. Effects of inhomogeneous brodening on one color pump-probe amplitude and phase profiles, (a) shows the ab- 
sorption lineshape function for the homogeneous case (dashed line), and with Gaussian inhomogeneities of aj = 300 cm~^ 
(solid line) and aj = 600 cm~^ (dash dot line). The amplitude and phase profiles for the ground state case are shown in 
panels (b) and (c), for the three lineshapes. The corresponding excited state amplitude and phase are shown in panels 
(d) and (e). 



In Fig. 6, we consider a homogeneous lineshape with Fg = 300 cm ^ and a single undamped 220 cm 
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mode with coupling S = 0.05. Pump-probe signals are calculated for T = 300 K using both the homogeneous 
lineshape and with Gaussian inhomogeneities included. The top panel in the Fig. 6 shows ^i{lu) with three 
values of the inhomogeneous width ct/ — 0, 300 cm^^, and 600 cm~^. The lower panels show the amplitude 
and phase of ground and excited state signals for the three cases assuming a pulse width of Tp = 50 fs. 

It is seen that while the amplitude dip and phase discontinuity of the ground state signal are dramatically 
removed by the inhomogeneity, they exhibit a red-shift for the excited state. The shift of the excited state signal 
arises from the asymmetry of the excited state amplitude: since the amplitude of the homogeneous stimulated 
emission signal on the red side is much smaller than on the blue side of Sly, an effective cancellation of the 
positive and negative lobes of the signed amplitude will occur towards the red. 




FIG. 7. (a) Equilibrium absorption spectrum of Mb at T = 300 K (open circle) shown with the theoretical fit (solid 
line) based on a prior model for the asymmetric deoxy Mb lineshape. The underlying homogeneous lineshape (before 
convolution with inhomogeneities) is also shown (dotted line), (b) Experimental Soret band phase profile for 220 cm~^ 
mode coherence from dcoxy Mb (filled circles) is shown with the theoretical fit based on an inhomogeneous model (solid 
line) and a homogeneous model (dashed line) using the asymmetric lineshape plotted in (a), (c) Experimental phase 
profile for the 220 cm~^ mode in MbNO for a range of carrier frequencies across the Soret band is shown with a theoretical 
fit based on simple model for reaction driven coherence in MbNO. 



1. Example: Inhomogeneous broadening in myoglobin 

The characteristic features of the ground and excited state amplitude and phase profiles can be used as 
indications of inhomogeneous broadening. As an example, consider the heme protein deoxy myoglobin (Mb), 
which is an oxygen storage protein found in muscle cells. Mb possesses a highly asymmetric and broad absorption 
spectrum (Soret band) in its ligand-free, high spin (deoxy, S=2) state. In previous studies, the Soret band of Mb 
was modeled using a non-Gaussian distribution of electronic energy levels, ascribed to disorder in the position of 
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the central iron atom of the porphyrin ring ]7l| , |75| ] . In pump-probe and resonance Raman experiments on Mb, 
the Soret excited state is found |7^j35|] to be very short hved (< 30 fs). Thus, the detected pump-probe signals 
for deoxy Mb are presumably due to ground state coherence driven by impulsive stimulated Raman scattering. 
Here, we focus attention on the Fe-His mode at 220 cm~^ which is a prominent mode in the resonance Raman 
experiments on Mb. 

In the top panel of Fig. 7, we plot the experimental absorption spectrum of Mb along with a theoretical 
fit based on the previously proposed model for the inhomogeneous distribution ||7l|,|7^. All resonance Raman 
active modes are included in the model, with coupling strengths determined from their absolute resonance 
Raman cross-sections |7^. The pump-probe calculations are carried out using the effective linear response 
approach using a realistic pulse width of 70 fs. For simplicity, only the 220 cto~^ mode {S = 0.05) was taken 
to be displaced from equilibrium. In Fig. 7(b), we plot the phase profile of the 220 cm~^ mode observed in open 
band (one color) pump-probe measurements on Mb across the absorption maximum. The theoretical prediction 
for the phase of the 220 cm~^ mode in the ground state is shown for both the inhomogeneous lineshape model 
as well as for a homogeneous model using the observed lineshape. 

It is clear from Fig. 7(b) that the theoretical prediction using the inhomogeneous lineshape is in good agree- 
ment with the experimental data. Note that the phase approaches the off-resonant limit of (7r/2) much more 
rapidly on the red side than on the blue side of the absorption maximum. This refiects the strong asymmetry 
in the absorption spectrum. The sharp phase jump between — (7r/2) to +(7r/2) predicted by the homogeneous 
model is clearly not observed in the experimental data. Also, The amplitude of the 220 cm^^ mode does 
not show a clear dip near absorption maximum, contrary to the prediction of the homogeneous model. These 
observations strongly support the presence of significant inhomogeneous broadening in deoxy Mb, as well as 
the ground state origin of the vibrational oscillations. If the 220 cm~^ coherence originated from the excited 
state rather than the ground state, an amplitude dip and phase fiip should be observed but red-shifted from 
the emission maximum as discussed above. 



D. Reaction driven coherence: 

Coherent oscillations in MbNO 



We now apply the effective linear response approach to coherence driven by a chemical reaction of the type 
shown in Fig. 1(b). The effective linear response function for a two level system consisting of the product ground 
and excited states |/) and |/') was discussed in section III.B.3. The nuclear density matrix for the state |/) is 
described by a displaced thermal density matrix. The displacements are the first moments of the vibrational 
coherence driven by a Landau-Zener surface crossing between the non-radiative states |e) and |/) presented in 
Appendix-B. The pump-probe signal arises from the vibrationally coherent state |/). As discussed earlier, we 
are concerned with the situation where the vibrational mode Q is not coupled to the .g — > e optical excitation, 
but is triggered solely by the non-radiative chemical reaction between |e) and |/) that follows pump excitation. 
For example, Q could represent the 220 cm~^ Fe-His mode in the reactive sample MbNO. The 220 cm~^ mode 
is not resonance Raman active in the reactant ground state {\g)) but is nevertheless observed in pump-probe 
photolysis signals. 

The spectral lineshape fun ction also plays an important role in the multidimensional Landau-Zener theory 
As specified in Eq. (B2) of Appendix-B, the transition rate between the states |e) and |/) is obtained by 
taking the zero frequency limit of the spectral function for the e ^ / transition i.e., = 0,QqI{R{s))). 

Here, r2QQ(i?(s)) is the time dependent energy gap between the |e) and |/) states as shown in Fig. 1(b) and 
R{s) is the time dependent semi-classical reaction coordinate. As the reaction coordinate sweeps through the 
crossing point, nll{R{s)) passes through zero. The reaction rate C/(i?(s)) is integrated as in Eq. ( pi] ) to obtain 
the time dependent quantum yield Pef{t). In the present discussion, we hold to the assumption of constant 
reaction velocity at the crossing point. This means that the energy gap I^qq (i?(s)) is decreasing at a constant 
rate /3. 

Our aim here is to illustrate the analysis of pump-probe signals associated with the reaction dynamics, and 
we restrict our attention to the role played by the reaction velocity /3. We refer the reader to reference pl| | 
for the details of the role played by the various reaction specific parameters in the population and coherence 
dynamics. In Fig. 8(a), we simulate the spectral function for the non-radiative states. The coupling of the 
220 cm^^ mode is taken as S* = 0.5 corresponding to the 0.1 A Fe-His bond length contraction associated with 
the loss of the NO ligand. An over damped low frequency bath mode is also included to broaden the spectral 
function. In the panels directly below the lineshape, the time dependent population transfer Pef{t) and the 
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first moment dynamics Qf{t) in tlie product state are plotted. Qf{t) is obtained botii by direct integration 
using Eq. (B3) and the initial first moments using Eqs. (B4-B6|). A reaction velocity of /3 = 2000 cm~^ /ps and 
P = 10000 cm~^ /ps are considered in Figs. 8(b-e). The initial vertical energy, corresponding to the separation 
of the quintet and singlet energy levels in MbNO, is taken to be nll{0) = 2000 cm . The connection between 
the time dependent population evolution and the first moment dynamics in the product state is clearly exposed 
in the simulations. We see that for smaller /?, the reaction proceeds in a smooth fashion, with a complete 
transition to the product occurring immediately after t — \ps. The slow passage through the reaction vertex 
causes the product state oscillations to be effectively diminished. For large /3, the product electronic state 
population changes more abruptly, and the initial amplitude of the oscillatory motion is considerably larger, 
approaching the equilibrium displacement of the reactant well with respect to the product (|Ae^| = 1). 

In Figs. 8(g-h), we simulate the amplitude and phase profiles of the pump-probe signals associated with the 
reaction driven oscillations in the electronic state |/), using the inhomogeneous equilibrium lineshape of Mb 
shown in Fig. 8(f). Note that only the product state electronic transition / ^ /' is inhomogeneosly distributed. 
The probe pulse detects the modulation of the / — > /' electronic transition by the reaction driven oscillations. 
Since the initial conditions of the wave-packet on the state | /) are independent of the pump wavelength, a sharp 
jump in the phase is observed accompanied by a dip in the amplitude at the absorption maximum nl. The direct 
correspondence between the amplitude of the vibrational oscillations and the amplitude of the signal is evident 
from the simulations. The slower reaction exhibits a much weaker signal amplitude than the faster reaction. 
The signal phase directly yields the initial phase (pf of the wave-packet according to Eq. (44b). Furthermore, 

if we let 7 = and Pef (t) = 5{t — tr) in Eqs. (B5) and (^), we have (pf = —ujoth, tr = rigQ(0)//3 being 
the reaction time. Thus, when /3 is increased, tr decreases and the initial phase continually changes through 
cycles of 27r. The phase approaches zero for an infinitely fast reaction. It is interesting to note from Fig. 8(c) 
and (e) that the effective initial momentum of the wave-packet (due to a non-zero (/S/) reverses its sign as /? is 
increased from 2000 cmT^ jps to 10000 cmT^ jps. Compare Fig. 3(c), where the momentum impulse has the 
same direction across resonance. 

As an example, we consider the Fe-His mode in the reactive sample MbNO in Fig. 7(c). In contrast to the 
non-reactive Mb, the phase of the 220 cmT^ oscillation in MbNO exhibits a sharp jump at a wavelength that is 
red shifted with respect to the Soret absorption maximum at ~ 23000 cmT^ (435 nrasoi the Mb product. Also, 
the amplitude of the 220 cmT^ oscillations vanishes at this detection wavelength js^]. These observations are 
consistent with the suggestion that the Fe-His mode is driven by the chemical reaction associated with the e f 
electronic transition as described in Fig. 8. The 4 nm red shift of the phase flip is probably associated with the 
initial red shift of the photo-product absorption band centered roughly at 439 nm (^ 22780 cm~^). The fact 
that the phase flip occurs precisely at this wavelength is consistent with the idea that the 220 cm^^ oscillations 
in MbNO are associated with the initial Mb photo-product that is created following photodissociation of the 
NO Hgand. 

It is also seen from Fig. 7(c) that the phase of the open band 220 cm~^ oscillations is near or tt. This 
suggests that the initial phase of the 220 cm~^ mode coherence is also close to 0. A zero initial wavepacket 
phase is realized if the reaction time tr is a small fraction of the 150 fs oscillatory period. The wave-packet 
is then placed almost instantly (^ 20 fs) on the product state potential surface. In the Landau-Zener theory, 
the reaction time depends both on the reaction velocity (3 and on the initial vertical energy gap ^IqI{R{0)), 
suggesting that either the gap is small or the velocity is large (or both) in the MbNO reaction. As discussed 
earlier, this corresponds to a population evolution that behaves like a step function. 

On the other hand, a zero initial phase of the wavepacket could also be an indication of direct pump pulse 
preparation of the product state coherence. For example, when riQQ(i?(0)) is very small, the pump pulse 
interaction may directly excite an adiabatic mixture of the |e) and the |/) electronic states. The coherence 
induced along the Q coordinate will then have a phase of zero (no initial momentum), since it originated from 
the excited state directly prepared by the pump pulse. This possibility is beyond the scope of the present 
treatment, which assumes that the pump pulse and chemical reaction can be separated. We will address this 
interesting question in future work. 

Finally, we point out that the absolute phase of the 220 cm~^ mode in MbNO indicates that the Fe-His 
bond in Mb (|/)) contracts upon optical excitation to Mb* (|/')), as depicted in Fig. 1(b). This is based on the 
knowledge that the Fe-His bond also contracts upon loss of the NO ligand; i.e. upon the non-radiative transition 
between |e) and |/). The wavepacket induced in |/) is therefore placed at positive displacements with respect 
to the equilibrium position. It is clear from Fig. 1(b) that if the Fe-His bond expands upon photo-excitation 
to I/'), exactly the opposite phase behavior from that observed in Fig. 7(c) would be predicted. This confirms 
earlier predictions [|78[ and also demonstrates that the phase of reaction driven pump-probe signals can provide 
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information on the sign of the optical electron nuclear coupling, provided the initial ground state equilibrium 
changes are known. 
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FIG. 8. Application of the effective linear response approach to reaction driven coherence, (a) Spectral lineshape 
function for the non-radiative e —f f electronic transition. Panels (b) and (d) show the time dependent population 
evolution of the product state for reaction velocities of P — 2000 cm~^ /ps and (3 — 10000 cm~^ /ps. Panels (c) and (e) 
show the corresponding me an p osition of the 220 cm~^ mode coupled to the reaction. The calculations using the effective 
initial conditions Eqs. (B4-B6) are shown in circles. The results of a direct integration as in Eq. (B3) are shown as a 
solid line. The non-radiative coupling strength was fixed at J = 120 cm~^ in both the calculations, (f) Spectral lineshape 
function for the product state (identical to the Mb lineshape in Fig. 7(a)). (g) and (h) show the pump-probe open band 
amplitude Ag (solid line) and phase <j)g (dashed line) profiles arising from the reaction driven 220 cm~^ coherence in the 
product state, for (3 — 2000 cm~^ /ps and (3 — 10000 cm~^ /ps. 
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E. Overtone signals 



1. First moment modulation (Coherent state) 

Up to this point, we have focused on the oscillatory pump-probe signal at the fundamental frequency of 
the optically coupled mode. The presence of higher overtones is expected even for harmonic modes, since the 
observed signal is not directly the mean nuclear position oscillating at the fundamental frequency. The pump- 
probe signal arises from a dynamic modulation of the equilibrium correlator as expressed in Eq. ( |3^ ) and jSq). 
Intuitively, if we picture a two electronic level system, the overtone signal can be attributed to the "double 
pass" of the wave-packet across the probing region The strength of the overtone signals due to the 

double pass is determined by the higher powers of the product AgA. This can been seen by expanding the 
exponential in the correlation function in Eq. (|3^). For weak field intensities, the pump induced ground state 
amplitudes are much smaller than the optical coupling strengths js^,^,^, which are often much less than 
unity in large multimode systems. Thus, the overtones generated by the first moment modulation of the pump 
induced ground state wave-packet are much weaker than the fundamental. 

On the other hand, when the coherence is induced by a non-radiative surface crossing or by intense pump fields 
|^ , |63| , the amplitude of the coherent motion can be comparable or larger than the optical coupling strengths, 
so that the strength of the overtones is expected to be larger. The effective linear response approach using 
displaced wave-packets allows a calculation of pump-probe signals for arbitrary values of the initial amplitude 
Ag. In order to estimate the strength of overtone signals as a function of Ag, we consider the model lineshape 
of Mb in Fig. 8(f). The overtone signals of a low frequency mode at 40 cm~^ coupled to the optical transition 
are studied in Fig. 9. The top panel of Fig. 9 shows the profiles of the fundamental and the first two overtones 
of the 40 cm~^ oscillation, assuming — = 1, and Sa = = 10. The contrasting behavior of the 

overtone and fundamental amplitude profiles is notable from the figure. The fundamental signal amplitude dips 
near the resonant maximum fly and peaks near the shoulders of the absorption spectrum. The amplitude of 
the first overtone peaks near Qy and vanishes close to the frequency where the fundamental dips. Thus, we can 
say that the amplitude profiles of the fundamental and the first overtone behave roughly as the first and second 
derivatives of the equilibrium absorption spectrum. Similarly, the second overtone at 120 cm~^ behaves as the 
third derivative of the absorption spectrum, vanishing at three different carrier frequencies across the absorption 
band. Such characteristics are useful for the assignment of overtone oscillations in pump-probe spectroscopy 

It is important to note the significant reduction in the amplitude of the overtone oscillations compared to the 
fundamental. To study the strength of overtone signals over a wider range of initial displacements, we plot in 
Figs. 9(b-d), the amplitude of the first three harmonics for three different carrier frequencies near resonance. 
The initial wave-packet amplitudes are varied over three orders of magnitude and Sa is fixed to be unity. It 
is seen that for detection at the wings of the absorption ujc — ^ci and Wcs, the amplitude of the overtones is 
more than an order of magnitude weaker than the fundamental. This is true even for amplitudes as large as 
Sa = 100 (or Ag K, 14). Near band center a;c2, the fundamental signal itself is very small. The first overtone 
is the dominating contribution even for very small displacements. The almost linear behavior of the overtone 
amplitudes on the log-log scale suggests a power law dependence of the overtone amplitudes on the wavepacket 
displacement. 

2. Second moment modulation (Squeezed states) 

While the contribution of the double pass of the wave-packet to the intensity of overtone signals is negligible 
for field driven coherence (see previous section) , higher moment modulations of the pump induced wave-packet 
p6| (squeezing) can contribute more significantly to overtone intensities. Overtones signals calculated using the 
full third order response approach are found to be much stronger than those calculated with the effective linear 
response approach using the first moment dynamics. This discrepancy is indicative of the role played by field 
induced squeezing in the ground state signals. Apart from the laser pulse interaction, squeezing of the nuclear 
wave-packet can also be induced by purely geometrical effects, such as a change in the curvature of the nuclear 
potentials that are involved in the interaction | |79[ i.e. quadratic coupling. 

We discuss squeezing in the context of chemical reactions. Introducing quadratic coupling between the 
potentials of the reactant and product surfaces |e) and |/) shown in Fig. 1, we calculate the probe response to 
displaced and squeezed dynamics on the product state potential well. For the initial displacement, we continue 



24 



to use the first moments of the Landau-Zener coherence presented in Appendix B. These were derived assuming 
linear coupling between the reactant and product surfaces. To obtain the initial conditions of the squeezed 
state, we adopt a phenomenological view. We assume that the effect of the reaction is to instantly place the 
thermal equilibrium distribution of the reactant on the product state well. The resulting non-stationary nuclear 
state is a displaced-thermal density matrix that is also squeezed due to the incommensurate thermal widths 
on the two surfaces. In Appendix-D, we have derived an analytic expression for the effective linear response 
correlation function that describes both the coherent and squeezed dynamics in the ground state. 

The central quantity governing the squeezed state dynamics is the ratio of the reactant and product state 
frequencies r = iuje/^f)- It is evident from the effective linear response correlation function Eq. (D6) that 
squeezed states modulate the equilibrium absorption correlator at even harmonics of the fundamental frequency. 
We therefore do not expect the fundamental pump-probe signal to be affected significantly by squeezed wave- 
packet dynamics. In Figs. 9(e-g), we simulate the first three harmonics of the 40 cm^^ oscillations, when the 
wave-packet is both displaced and squeezed (The amplitude and phase profiles in the presence of squeezing 
behave similar to those in Fig. 9(a) and are not plotted). In the simulations, the initial displacement of the 
wave-packet on the |/) state and the optical coupling are fixed at Sa — Sa — 1. The squeezing ratio r is varied 
over a range of 1 to 8. It can be seen that fundamental oscillations are practically unaffected by variations 
in r. The dramatic effect of squeezing on the overtone amplitude is evident, particularly so for the detection 
frequency U!c2 near the absorption maximum. Here, for squeezing ratios larger than 6, the overtone signal is 
comparable to the corresponding fundamental signal at the wings of the absorption band (wd, Wcs)- 
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FIG. 9. Simulations of overtone signals from the ground state of Mb using displaced and squeezed states, (a) 
Amplitude profiles of the fundamental at 40 cm~^ (solid line) and the first two overtones at 80 an~^ (dashed line) 
and 120 cm~^ (dotted fine), assuming an initial displacement of Sa = = 10, and optical coupling strength of 

S'a = = 1, for the 40 cmT^ mode. Panels (b-d) show the amplitudes of the first three harmonics over a wider 

range of initial displacements, for three carrier frequencies (lo^\,u)ci^i^cz) across the absorption maximum. The corre- 
sponding amplitudes of the three harmonics for displaced and squeezed states is plotted in panels (e-g) , for a wide range 
of frequency ratios r = uOe/ujf between the reactant (|e)) and product (|/)) electronic states. 
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The simulations presented in Fig. 9 provide an overall picture of overtone intensities in pump-probe spec- 
troscopy, within a harmonic model with coherent and squeezed initial states. It can be concluded that the 
fundamental signal is almost always the dominant contribution for reasonable values of wave-packet displace- 
ment and squeezing. However, when quadratic coupling is also introduced between the electronic levels involved 
in the probe interaction, i.e. when ujf ^ LOfi^ the overtone signals are expected to be significantly enhanced. 
In this case, the strength of the overtone oscillations induced by the curvature change (between the ground 
and excited state potentials) is coupled to the amplitude of the wave-packet motion; larger amplitudes imply 
that steeper differences in curvature are probed. In contrast, the overtone intensities generated by squeezing 
the initial wave-packet are dependent only on the frequency ratio. A rigorous treatment of the effective linear 
response functions for a quadratically coupled two level system will be presented elsewhere. 



V. SUMMARY AND CONCLUSIONS 



In this work, we have presented a study of the amplitude and phase excitation profiles in femtosecond 
coherence spectroscopy using an effective linear response approach. In combination with a rigorous moment 
analysis of the doorway state 1 46 1 , the linear response formalism is found to be in excellent agreement with the 
more complete third order response formalism [ p8[ . Furthermore, the representation of non-stationary states as 
displaced thermal states provides a general framework to analyze vibrational coherence induced by non-radiative 
transitions. We have demonstrated this approach by using the first moments of a quantum coordinate coupled 
to a Landau-Zener surface crossing. 

The effective linear response functions are analyticaband yet retain the correlation function based description 
(in contrast to the vibronic eigenstates approach implicit in the third order response formalism. The 

analytic expressions offer a significant reduction in computation times. Calculations of multiple carrier frequency 
excitation profiles for arbitrary pulse widths and temperatures are made possible, without the need for semi- 
classical approximations. A separate moment analysis is also useful in providing clear physical insight into 
pump-probe signals. Ground state signals generated by impulsive Raman processes exhibit contrasting phase 
and amplitude behaviour from excited state and reaction driven signals. 

The ability to incorporate experimentally measured absorption lincshapes (and the Kramers-Kronig deter- 
mined dispersion lineshapes) into the pump-probe calculations is another important aspect of the present work. 
A knowledge of pump-probe excitation profiles and the measured absorption cross-sections can be used to ex- 
tract parameters relevant to a specific mode that is active in pump-probe signals. The use of experimental 
lineshapes eliminates the need for modeling the multimode vibronic mechanisms that cause line broadening. 
This is analogous to "transform" techniques in resonance Raman spectroscopy [p8|- |70|j5^ . Transform techniques 
are a powerful means for extracting mode specific information using the measured Raman excitation profiles 
and the optical absorption profiles on an absolute (cross-section) scale. 

A study of inhomogeneous broadening in pump-probe signals shows that the ground state amplitude and 
phase profiles are smeared out by the inhomogeneous broadening. The excited state amplitude dip and phase 
flip are red shifted with respect to the peak of the emission lineshape. A pump-probe analysis of the phase 
profiles of the 220 cm~^ Fe-His mode in Mb across the Soret band was presented, and was based on a previous 
model for the strongly asymmetric Soret lineshape . The analysis indicates both the ground state origin of 
the oscillations and the strong inhomogeneous broadening present in Mb. In contrast to the nearly constant 
phase profile in Mb, the 220 cm~^ mode phase in the reactive sample MbNO shows a clear tt phase change near 
the peak of the transient photo-product. These results suggest that the 220 cmT^ mode coherence in MbNO 
arises from the photo-product Mb. The mode is presumably triggered by the rapid curve crossing occurring 
during the MbNO photolysis reaction |Q. The absolute phase of the oscillations suggest that the chemical 
reaction occurs on a very fast time scale (^ 20 /s), with the possibility of direct pump pulse involvement at 
the crossing region in the preparation of the product state coherence. 

Although we have not discussed vibrational relaxation in detail, the non-stationary linear response correlation 
functions in Eqs. (jsj) and (^5|) can readily incorporate damped wave-packet motion. The initial first moments 
for the damped oscillator case must however be used instead of the undamped expressions in Appendix-A. In 
the analysis of the oscillatory amplitude and phase profiles, which has been the main subject of the present 
work, the neglect of damping is justified because the observed modes are weakly damped. However, both under 
damped and over damped modes contribute to the non-oscillatory background (offset) in the pump-probe signals 
[ 44 , ^j20| . These offsets can potentially carry important information concerning the equilibrium potential shifts 
[20|. It is immediately clear from Eqs. (^4|) and ( ^5|) that the modulation of the equilibrium correlator by 
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strongly over damped or diffusive motion manifests itself as a dynamic i.e. shifting absorption lineshape in the 
ground state response . Overdamped motion is manifested in the excited state response through phenomena 
such as fluorescence Stokes-shift dynamics . A rigorous analysis of these processes can be readily carried 

out using the effective linear response approach formalism presented here. 

We have also discussed overtone signals in pump-probe spectroscopy using both displaced and squeezed 
initial states. The contrasting amplitude and phase behavior of the various harmonics are potentially useful for 
overtone assignments. It is found that for reasonable values of squeezing ratios and wave-packet displacements, 
the fundamental signal is the dominant contribution. 

Finally, we mention that the present development can be readily extended beyond the displaced oscillator 
model. Effective linear response functions can be readily calculated for the case of quadratic electron-nuclear 
coupling (between the radiatively coupled levels) using standard techniques of quantum field theory. Non- 
Condon effects can also be incorporated 
These calculations will be detailed elsewhere. 
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into both the pump (moment analysis) and the probe steps. 
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APPENDIX A: FIRST MOMENTS OF PUMP INDUCED VIBRATIONAL COHERENCE 



Consider the ground state doorway 5'pg in Eq. (p^). Let the pump field be of the form Ea{t) = EQG{t) cos{LUct) 
where G{t) is a dimensionless envelope function, and Eq is the electric field strength of the pump pulse. General 
expressions have been derived for all the moments of Q and P for the density matrix Spg using moment generating 
functions, and related to equilibrium lineshape functions [^^. The first moments of Q and P are found to be 
(in dimensionless units): 



/•oo 

/ dwGp(w - Wc, -wo)A$/(w), (Al) 

"'0 



\f,g,\'EU2n+l)A 
" Snh'N'g 

P90 - r' 2 ° / dioGpiio ~ uj,, -c^o)A$fl(c^). (A2) 
87r?i N'g Jo 

Here, f^'g is the net ground state population after the pump interaction and <l>/(w) and $^(0;) are the imaginary 
and the real parts of the complex equilibrium lineshape function defined in Eq.(|2^). Gp{uj — Wc, — wq) is the 
product spectral function defined in Eq. (^TJ). The action of the operator A (not to be confused with the 
dimensionless excited state potential shift A) is to generate differences: A$/(a;) = ^i{io) — $/(w — ojq). With 



the initial position and momentum increments given in Eqs. (Al) and (|A3), the time dependent first moment 



dynamics in the ground state is given rigorously as the solution to the damped oscillator equation: 

Qg{s) = Age-''\'\ cob{u,s + (^g), (A3) 



where 7 and u^, = \/ are the damping constant and effective frequency of the mode. The amplitude 



and phase are given hy Ag ^ ^ Q'^g^ + Pg^ and <Pg = - tan ^ [Pgo/Qgo]- 

The various moments of Q and P for the excited state doorway (5pe in Eq. (^6|) can also be obtained. It is 
found |^6|] that the excited nuclei receive no initial momentum i.e., Peo = 0. The mean position (with respect 
to the excited state potential equilibrium A) is given by: 

QeO = , ^2.r, / dujGp{u) - LJc,~UJo) 

Ann A*; Jo 

^/(w - Wo) -nA$/(w) , (A4) 

where N'^ is the electronic population in the excited state {Ng + = 1). Thus, the amplitude for the excited 
state coherent motion is given by A,, — \Qeo\, and the time dependent first moment is 

Q^{s) = Aee-^l'^l cos(a;^s). (A5) 



APPENDIX B: FIRST MOMENTS OF COHERENCE DRIVEN BY LANDAU-ZENER SURFACE 

CROSSING 

We consider a product electronic state |/) that is coupled non-radiatively to the reactant excited state |e) 
as in Fig. 1(b). The Hamiltonian of the problem consists of the 2x2 lower block diagonal matrix of Hjsm in 
Eq. (^6|). The system is assumed to be initially (t = 0) in thermal equilibrium along the quantum mechanical 
degrees of freedom of the state |e). The non-radiative couphng J will induce a transition to the state |/). 
The time-dependent surface crossing along the classical reaction coordinate R will be accompanied by quantum 
mechanical tunneling between the vibrational levels of He and Hf. The time dependent quantum yield for the 
transition e ~* f can be shown to be: 

Pef{t) = 1 - exp | - ^ dsg{R{s))^ . (Bl) 
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Here, G{R{s)) is the Fermi Golden rule transition rate for making an electronic transition from the state |e) 
the electronic state |/) and is expressed as |2^: 



to 



5iRis)) 



(B2) 



illl{R{s)) = Uf{R{s)) — Ue{R{s)) is the time-dependent energy gap between the e and / states, (uj) is the 



/ electronic transition, obtained using Eqs. 




with 



spectral absorption lineshape function for the e 

^00 replaced by Qll{R{s)). In the single mode hmit, i.e. no Q coupling, Eq. ( |B2D yields a delta function. The 
time dependent quantum yield Pef{i) will then be a simple step function in time. The presence of additional 
degrees of freedom smears out the curve crossing process and can be completely described by the zero frequency 
component of the spectral lineshape function <^'^/ {lj). 

The time dependent surface crossing along the R coordinate induces oscillations along Q, which is initially 
in thermal equilibrium. It can be shown generally that the mean position of the oscillator in the product state 
obeys 



Qfit) 



dse '^'^ sm{ujys)Pef{t — s). 



(B3) 



Here, 7 is the vibrational damping, Ae/ is the dimensionless shift between the e and / oscillators, and ojy = 
a/o^o ~ 7^- The above equation rigorously describes the first moment of the nuclear dynamics on the product 
state well, with the one qualification that it is valid only at times t when the system has made its transition 
into the product state, i.e. Pef{t) = 1. When Pef{t) < 1, Q f{t) consists of oscillations in both the reactant and 
product surfaces and is more difficult to interpret. We can recast Eq. (B3) to a form that readily yields the 
initial conditions for the probe i.e. detection step. Integrating Eq. (B2) by parts, we get in the long time limit 
(assuming Pe/(oo) = 1) 



where 



Qj:{t) = Aef + Afe '^*cos(w„t + v'/), 

/OC 
-00 

Afsm{ipf) = Ae/ / dt' Pef (t')e'''' [sin(w.„f') + {-//iVy) cos{ujJ')] . 



(B4) 

(B5) 
(B6) 



Eqs. (B5) and (B6) determine the amplitude and phase of the oscillatory motion induced by the surface crossing. 
They can be readily incorporated into the probe step using a displaced thermal representation. The above 
expressions depend on the reaction parameters (J, /3, figg (i?), etc.,) through Pef{t) and Eq. (Bl). Extension of 
the multi-dimensional Landau-Zener theory to incorporate non-stationary initial conditions along the Q degree 
of freedom can be carried out using displaced thermal density matrices, with the displacements obtained from 
the excited state first moment presented in Appendix A. 



APPENDIX C: NON-LINEAR RESPONSE FUNCTIONS 

When the ground state doorway 6pg in Eq. ( p3| ) is substituted in the expression for Cg{t,t^) in Eq. (p^), 
we are left with the thermal average of a four-time correlation function. The average can be evaluated exactly 
using a second order cumulant expansion and the result is 

2 poo .t2 

Cg{tM)^-^ dt2 dtiEa{tl)Ea{t2)[Ri{tMM,tl)+R3{t.h,t2,tl)]. (CI) 

1^ J —00 J —OG 

Here, R3 and R4 are the non- linear response functions given by ||83| , ^ , p8| , 

Rsit, h,t2,h) = Hit - h)n{t2 - ti)e-*"™(*-*^-*^+*i)e-9(*-*^')-3'(*=^-*i)+-^^(*'*i'*^'*^\ (C2a) 
-R4(i,ti,i2,i3) =H(t-t3)H(i2 -ii)e-'"™(*-*^+*'-*^)e-9(*-*3)-fffe-ti)-/2(t,ti,t2.t3)^ (C2b) 
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where we have defined 



f2{t, hM.h) = g{t -h)-g{t- ts) + g{t3 ~ tz) - gits - h). (C3) 

The function 7i(t) accounts for homogeneous broadening and is typically of the form e^'^^'*', being the 
homogeneous damping constant. In these expressions, t2 — ti reflects the time interval when the system is in an 
electronic coherence during the pump field interaction, whereas t — t^ is the time interval when the system is in 
an electronic coherence during the probe interaction |2^. Thus, Hit) appears separately for the pump and the 
probe interactions. 

Turning to the excited state response, substitution of 6pe in Eq. ( p6| ) into the expression for Ce in Eq. ( ^5| ) 
and using a second order cumulant expansion, we find 

2 j-oo .t2 

Ceit,t3)^^ dt2 dhEa{ti)Ea{t2)[Rl{tMMM)+R*2{tM,t2M)]- (C4) 

Here, Ri and i?2 are the non-linear response functions given by ||26| , ^8|| 

i?i(i,ii,t2,t3) ==H(i-t3)H(i2 -ti)e-*"'"'(*~*-'+*'"*^'e-f'(*-*^')-9(*^-*i)--'^i(*'*i'*^'*3', (C5a) 
R2{t,ti,t2,t3) t3)H(t2 -ii)e"'"™'*"*^~*"+*^^e-s'(*-*3)-s*(t2-ti)+/r(t,ti,t2,t3)^ (C5b) 

and we have defined 

hitMMM) = g{t - h) -g*{t~ h) " g{h - h) + g*{h - h). (C6) 

In the above expressions gis) is the harmonic oscillator correlation function [Eq. (^)]. 

Equations ( |Cl| ) and ( |C^ ) express the two time effective linear response correlation functions Cg^e{t-,tz) as 
as convolution of the well known non- linear response functions [||j2^ Rj (j = 1..4) with the two pump field 
interactions. The convolution with the pump laser fields implicitly takes into account the complete information 
about the non-stationary state of the system expressed by the density matrices Spg and dpe- It is worth noting 
here that both Cgit^t^) and Ceitjt^) can be factored into equilibrium and non-equilibrium parts as follows: 

Cg{t,t:i)^Kg{t-ts)C'{t,t3); Ce{tM) = K,{t - h)C',{tM). (C7a) 



where Kg{t — t^) and Kf.(t — t^) are given by Eq. (|27|). 

Inhomogeneous broadening due to a static distribution JF/(^^oo) can be incorporated into the above ex- 
pressions. The correlation functions Cg and Ce must be averaged with respect to the inhomogeneous dis- 
tribution. The electronic energy gap appears in the response functions Rj through the oscillatory factors 
g-inoo[(t-t3)±(t2-ti)] ^ ^]-^g averaging over Jv(i^oo) thus involves a simple Fourier transform and can be carried 
out independently of the time integrations. The resulting inhomogeneous response functions i?^^"* can then be 
expressed in terms of the homogeneous ones Rj in the following factorized form p8[ : 



R^l\{t, ti,t2, tj) — Ri^4:{t, ti,t2, t3)Ti{t — t3 + t2 — ti), (C8a) 
R^'lit, h,t2, h) = i?2,3(i, tl,t2, h)Ti{t -h-t2+ tl), (C8b) 



where J'lit) is the Fourier transform of the inhomogeneous distribution J-l{^l, 



APPENDIX D: EFFECTIVE LINEAR RESPONSE FUNCTION FOR A TWO LEVEL SYSTEM 

WITH SQUEEZED VIBRATIONAL STATES 

Consider two non-radiatively coupled reactant |e) and product |/) electronic states. Let the vibrational 
potentials be harmonic with linear and quadratic coupling, with frequencies Uf. and w/. Since two different 
vibrational frequencies are involved, we will first use dimensional units for clarity and will revert to dimensionless 
units later. Just before the reaction, the system is in thermal equilibrium on the electronic state |e) with the 

coordinate and momentum uncertainties {a\ ^IA^\ — IA\ ) given by 
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{Ti/2muje) (2nT 



i);4e 



[mujJi/2) (2nT + 1), 



(Dl) 



where Ht is the mean occupation number at temperature T. The position and momentum uncertainties in 
the product state are obtained from the above expression by simply replacing We by w/. If we picture the 
chemical reaction as instantly creating distributions of the above widths in the product state well, the initial 
non-equilibrium state in the product well will appear to be squeezed in the position and momentum coordinates. 
The non-stationary state can be generated from the thermal equilibrium density matrix of the product well by 
the squeezing operator defined as 



S{P) = exp 



;/3 



■In 



(D2) 



The present choice of the squeezing operator maintains the minimum uncertainty product. Arbitrary changes 
in the uncertainty product can be incorporated ||59t] by a complex value for the squeezing parameter /3, but 
is not considered here for simplicity. It is easy to verify that the unitary operator S{P) converts the thermal 
density matrix for the state |/), i.e. to the thermal density matrix for the state |e), i.e. The position 
and momentum uncertainties are also changed to the thermal values in the reactant state given in Eq. (Dl). 

Since the wave packet induced in the state |/) is in general both displaced and squeezed, we may represent 
the corresponding density matrix as 



With the above representation, we evaluate the correlation function similar to Eq. (|2 

1/2 

dimensionless units where the scaling factor is (h/mujf) ) 



(D3) 

We find (reverting to 



(D4) 



where we have used the invariance of trace under cyclic permutation to bring the unitary operators inside the 
time ordered exponential. Using the properties of coherent and squeezed states k^\ , we have 



Af cos(w/s -I- iff) + (aA(s) +a^A*{s)) 



(D5) 



where A(s) — (cosh(/3)e~*'^J"' + sinh(/3)e"^^^). Substituting this in Eq. (D4), we can remove the c-number 
displacement outside the time ordering. The resulting thermal average can be evaluated exactly using Wick's 
theorem. We finally obtain 



Cf{t,h) = Kf(t - t3)exp [{lAfA {s\n{ujft + ipf) - sin(w/t3 + ipf))] 
X exp [—R (2 cos{ujf{t + t^)) — cos{2ujft) — cos(2a;/i3))] , 



(D6) 



where we have defined 



A2 

R = — (2nT + 1) cosh(/3) sinh(/3) 



{2nT + 1)A2 (r2 - l) 



8r 



(D7) 



Here, r ~ uje/cuf and A is the dimensionless shift between the product ground and excited potentials. The 



equilibrium correlator is given by Kf{t — t^) 
as compared to Eq. (p8|): 



-in, 



ooe s(* where g{s) now takes a slightly different form 



A^ 

git - tg) = ^ [(2n^ + 1) (1 - cos{u;f{t - 13))) isin(w/(t - ta))] , 



(D8) 



where 



(2n^ + 1) = coth (hujf/2kBT) (cosh2(/3) + sinh^ (/?)) 



(D9) 
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Eq. (D6) expresses the non-equilibrium response as the modulation of the equilibrium correlator Kf{t — ts) by 
the coherent and squeezed dynamics of the nuclear motion on the ground electronic state |/) of the product. In 
addition to the exponential factor that describes coherent dynamics, squeezing introduces an additional factor 
that modulates at the even harmonics of the fundamental oscillator frequen cy tu f. It is seen that the strength 
of the overtone modulations arc determined by the factor R defined in Eq. (D7). We also note that the effect 
of a sque ezed initial state on the equilibrium part of the correlator is to change the effective temperature as in 
Eq. (pgh. 
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